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I. SUMMARY 


Theoretical and experimental studies have demonstrated the potential for 
aircraft capable of operating at hypersonic speeds within the atmosphere with 
an air breathing engine employing supersonic combustion. Duplication of true 
flight conditions in ground testing these propulsive systems requires a signi- 
ficant extension of existing facility capabilities and operating ranges. How- 
ever, the nature of the flow at these extreme conditions does not completely 
si mu late the actual flight conditions experienced by the Supersonic Combustion 
Ramjet (SCRAMJET) engine because the flowing gas stream Is not In complete 
thermal and chemical equilibrium. 

This research program Is a theoretical and experimental Investigation of 
the effect of nonequilibrium conditions upon the performance of combustors 
employing supersonic flows. Calculations and experiments are made regarding 
the effects on the Ignition of hydrogen of the nonequlllbrlum species (free 
radicals, atoms, water vapor, etc.) obtained using vitiated air. 

Results of this Investigation show that the nonequlllbrlum free-radical 
content from a supersonic vitiated air source will cause early Ignition of the 
hydrogen. An analysis of heated air expended from a high temperature source 
to test section conditions also Indicates that there Is sufficient free radi- 
cal content In the Incoming flow to cause early Ignition. Water vapor, an 
Inherent contamlnent In the generation of vitiated air, was found to reduce the 
ignition delay period under the experimental conditions considered. 
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II 0 INTRODUCTION 


Propulsion systems utilizing supersede combustion have recently been 
recognized as having strong potential for vehicles operating in the hypersonic 
flight regime. As flight velocities are Increased above Mach Stebet* '6 B ‘ the ; 
stagnation temperatures and pressures required for meaningful ground test • 
duplication becomes Increasingly difficult t© achieve, particularly, where large 
test sections and long experimental run tiroes are required, A typical super- 
sonic combustion ramjet engine (SCRAMJET) cruising at a Mach number of 6 at an 
altitude of 80,000 feet would encounter stagnation temperatures and pressures 
In the realm of 2500®R and 750 psla with air mass flows of 150 pounds per 
second. Dr. Antonio Ferri (1)* pointed out that the most urgent problem to be 
solved before a supersonic combustion engine could be perfected was the devel- 
opment of ground facilities for testing these engines. Within the present 
state-of-the-art of ground test facilities. It Is just about Impossible to test 
a complete engine under conditions duplicating all the parameters corresponding to 
to high Mach number flight, particularly for a time duration sufficient to 
investigate the time-dependent facets ©f the combustion process. In the super- 
sonic combustion problem as applied to the SCRAMJET, the chemical composition 
of the flowing medium, the mixing of the streams of fuel and air, the chemical 
reactions taking place, and the reaction rates involved are the Important fac- 
tors to be considered (2). The primary parameters requiring actual duplica- 
tion In the test cell are the static temperature, static pressure, composition 
of the gases, residence time In the combustor , and the velocity relationship 
between the mixing streams. Duplication of the severe environments associated 
with lea# altitude high Mach number flight regimes using state-of-the-art ground 
test facilities, is possible at Mach numbers up to approximately 7, but It Is 
difficult due to the size and operational cost of a clean air heater capable 
of delivering the necessary high temperature, pressure and quantity air mass 
flows. An alternate approach to this problem Is to generate a fluid medium 
labelled “vitiated air", a synthetic mixture of gases having the same percent- 
age oxygen content and molecular weight as true air. The most significant 
difference between vitiated air and true air Is the contaminating constituents 
such as water vapor and carbon dioxide emanating from the chemical reaction of 
generation. These contaminating constituents are responsible for uncertainties 
In the results, as compared to true air, when employing vitiated air for high 
temperature combustion research. On the ether hand, facilities employing 
vitiated air are able to operate at temperatures up t© approximately 4500 ®F 
and pressures up to 2-3000 psi together with high flow rate capabilities. 

This research program Is a theoretical and experimental investigation of 
the effects of nonequilibrium conditions — thermal and chemical — upon the 
performance of combustors employing supersonic flows. Specifically , the 
effects of free radicals, atoms and third-body species such as water vapor (a 
contaminating constituent In vitiated air) , upon the Ignition process are de- 
termined. Calculations and experiments are made regarding the effects on the 
Ignition of heated hydrogen of the nonequilibrium species obtained using 
vitiated air. 


*Numb@r in bracket indicates references In the Sibil ography. 



A. SCOPE OF PRESENT INVESTIGATION 

In the supersonic burning of gaseous hydrogen with air, there Is a finite 
distance required for the mixing and subsequent Ignition and combustion of the 
two gas streams. It Is a welt documented fact that the reaction rate of a 
chemical reaction Is directly proportional to the temperature and pressure, 
with the effect of temperature the more prominent (1). Thus, with a high ini- 
tial static temperature (T > 1200 # K) , reaction times are almost negligible and 
the mixing phenomenon of the two streams becomes the rate controlling mechan- 
isms. At low static temperatures corresponding to low flight Mach numbers , the 
chemical reaction times are slow and the reaction rate determines the combustor 
length. The present Investigation Is In the latter regime where the chemical 
reaction times are assumed significantly longer than the corresponding mixing 
times. 

In the ground testing of a SCR/WJET engine combustor, the working fluid 
must be expanded from an Initially high stagnation temperature regime (T & 
2200°K) to the relatively low temperature regime at the combustor entrance. 

Local conditions at this station are usually a fairly high velocity (6000 
ft/sec) stream at this relatively low static temperature (T » 1000 # K). During 
expansion' process free radical concentrations formed at the upstream stagnation 
conditions are frozen due to the fact that the gas residence times are much 
shorter than the reaction times at local temperatures throughout the expansion 
nozzle. This phenomenon known as "recombination freezing" results In nonequi- 
librium free-radical concentrations In the Incoming stream at the combustor 
Inlet. Initial free-radical concentrations In the incoming free stream have 
been found In previous analytical Investigations to have an Influenclal effect 
on the Ignition delay distance (3). 

In the formulation of the analytical portion of this Investigation, 
papers of several authors who investigated the Ignition delay mechanism of the 
hydrogen-air reaction systems were reviewed. The results which were felt to be 
germaln to the unique vitiated air-hydrogen reaction were incorporated Into the 
analytical model. Of particular concern were the effects of water vapor (an 
inherent constituent in vitiated air) , a substantial initial concentration of 
free radicals and relatively low static temperatures, on the combustion of 
heated gaseous hydrogen In a supersonic flowing stream. R. P. Rhodes (3) found 
that the initial free-radical concentration of the 1 neoml ng streain reduced the 
Ignition delay time by a factor of two for an order of magnitude Increase In 
Initial free-radical content. Vitiated air, Inherently In its generation has a 
substantial concentration of free radicals as will be shown later. Snyder et al 
(4) found that H 2 0, reduced the Ignition delay time an order of magnitude below 
that for dry air. Vitiated air has an appreciable amount of water vapor, 

(« 19% by weight). 

R. S. Brokaw (5) considered the Importance of Including H0 2 and H 2 0 2 reac- 
tions In low temperature ignition delay studies. These species , which were not 
Included In any of the other Investigations, were found analytically to Increase 
the Ignition delay In the temperature regime of Interest In this Investigation. 
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ess is assumed to occur in two steps 
es schematically the Ignition model, 
s are exponentially produced by: (1) 
tlal chain Initiating steps* and (2) the chain branching mechanisms of the 
blmolecular reactions. Throughout this period the temperature remains essen- 



period. The Ignition model assumption of constant temperature along with the 
simplification of constant pressure allows one to neglect the fluid-dynamic 
equations ami to solve only the differential equations describing the reaction 
kinetics. 


Most of the theoretical and experimental work concerning the combustion of 
hydrogen-air* is based on the assumption that complete thermal and chemical 
equilibrium exists at the combustor Inlet. In the shock tube work of Schott and 
Kinsey (6) and the standing detonation wave work of Nichols (7) , the combined 
streams are assumed In equilibrium Immediately downstream of the shock front. 



program for generating the specie concentration for the vitiated air. This pro- 
gram was developed by TRW Systems for the NASA (MSC). The program solves the 



the equations governing the rate of growth of free radicals and to solve these 
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FIGURE 1 SCHEMATIC OF THE IGNITION MODEL 







for species concentrations as a function of time using standard matrix tech- 
niques. The computer program used was a modified program supplied by the NASA 
Lewis Research Center. It is a simplification of a more general program devel- 
oped by the General Applied Science Laboratory (GASL). The modification in- 
cluded a constant ignition temperature simplification and the inclusion of 
several hydrogen-oxygen reactions (H0 2 , H 2 0 2 previously mentioned) which were 

analytically found to be important to the ignition kinetics at temperature 
below 1200°K. Inclusion of these species in the Induction analysis should re- 
sult in a more favorable comparison of experimental and theoretical ignition 
delays at low ignition temperatures. This computer program is labelled the 
Simplified Ignition Lag (SIL) program. 

Figure 3 is a schematic diagram for the analytical mechanism of the inves- 
tigation. Fifteen reactions are included in the SIL program. The main limita- 
tion of this program is that It falls to Include any NO or N0 2 reactions which 
were experimentally found to catalytically reduce the ignition delay time. 

Since the mechanism which accelerates the reaction is not understood at this 
time. It was not possible to include any NO or N0 2 reactions in a manner that 
would be meaningful. Results from the SIL program are species concentrations 
as a function of time. The ignition delay time Is determined from the semi- 
log plot of the hydroxyl radical concentration against time. This delay period 
is characterized by the region of essentially constant exponential growth of 
free-radical concentration. The end of the Ignition delay period is considered 
to be the time when the hydroxyl concentration deviates from constant exponen- 
tial growth as shown in the Figure 4. This definition for the end of the delay 
period is recommended by Hersch (9) who investigated the effect on ignition 
delay times for di fferent methods of identifying the termination of the delay 
period. 

Before applying the SIL program to the analysis of ignition delay with 
vitiated air , it was first necessary to establish the validity of the program's 
results. Since most of the theoretical and analytical work on Ignition delay • 
is concerned with real air and hydrogen, the verification was made for these 
gases. Initially a base case was selected for making this comparison and is 
shown with the base conditions also In Figure 4. As a starting condition both 
streams are assumed to be in chemical equilibrium at the assigned temperature 
and pressure prior to mixing. The combined stream is no longer in chemical 
equilibrium. In the bottom right hand corner are delay calculations of three 
other investigators for the base case conditions. In the next three figures 
are shown the effects of temperature, equivalence ratio and water vapor on the 
delay time under base case conditions. 

The effect of temperature on Induction time is shewn in Figure 5, in 
which hydroxyl radical concentrations as a function of time are presented at 
three separate temperature levels. The results show that a 100°K reduction in 
temperature over the base case just about doubles the ignition lag time. How- 
ever, a further reduction of the same magnitude, results in a factor of six 
difference in lag time over the base case. 

Figure 6, illustrates the effect of equivalence ratio on the ignition lag 
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FIGURE 3 ANALYTICAL MECHANISM FOR INVESTIGATION 




OH SPECIES CONCENTRATION, GRAM MOLES/CM 3 


(1) HYDROGEN-AIR 

(2) T » 1200°K, P * 1.0 ATM., ♦ « 1.0 

(3) EQUILIBRIUM INITIAL CONDITIONS 


x 10 - 39 pSEC 


APPROXIMATE DELAY CALCULATIONS 

SOURCE t j D 

(1) FERRI CG 24 m SEC 

(2) NICHOLES [6] 34 mSEC 

(3) RHODES QJ| 75 mSEC 


TIME, MICROSECONDS 


FIGURE 4 IGNITION DELAY FOR BASE CASE SPECIES 









FIGURE 6 EFFECT OF EQUIVALENCE RATIO ON DELAY TIMES 


for the base case species concentration. Three values for the equivalence ratio 
are presented on this figure. The results show that equivalence ratio has a 
very small effect on Induction times In the rang® (0.2 < § > 1.0). Momtchiloff 
(10) obtained similar results In a theoretical and experimental comparison of 
the effects of equivalence ratio on delay times. The analytical model of the 
Ignition process assumes instantaneous mixing for the two streams. This 
assumption Is justified on the basis that local mixture ratio variations during 
the actual mixing process have essentially no affect on Induction time. 

Production of vitiated air usually implies the presence of water vapor as 
a contaminant in the gas stream. The effect of water vapor on the Ignition lag 
time Is illustrated In Figure 7, which presents OH concentrations with time for 
the base case species and for a 10 percent water vapor additive. The 10 per- 
cent replaces an equal volume of air. Temperature, pressure, and equivalence 
ratio are fixed in this comparison. Water vapor has previously been shown to 
reduce the ignition time by a small amount as substantiated by Snyder et.al. (4) 

Several authors have derived either analytical or empirical relations for 
predicting the Ignition delay times for the hydrogen-air system. A comparison 
of their results with those obtained from the SIL program are presented in 
Figure 8. The correlations shown here are terminated at their respective tem- 
perature range of application. Note, that most of the correlations presented 
are plotted as straight lines. However, the SIL program results deviate consid- 
erably from straight line behavior especially at lew temperatures. The reason 
is that the SIL program Includes the effect of the H0 2 specie on the ignition 
times. This specie has often been observed In hydrogen-oxygen flames as an 
intermediate product that disappears rapidly as the reaction progresses. 

Brokaw (5) Is the only other Investigator to Include this specie and the point 
shewn includes his postulated mechanism by which the H0 2 specie influences the 
delay time. 

A comparison of the SIL program with experimental data for Ignition delay 
at low temperatures is shown In Figure 9. Above 1000°K the SIL program pre- 
dict delays close to those of Momtchlloff (10) and Schott (6) who Is Indicated 
on the graph by the square symbols. The SIL also tends to follow the basic 
trends of Snyder's shock tube data (4) which Is identified by the circles and 
triangles at the lower temperatures. All In all, the general agreement of the 
SIL program with experimental data tends to support the argument that Ignition 
delays Increase significantly at low temperatures due to the self-inhibition 
of the reactions by H0 2 formations. 

Examination of the ignition model en^loyed In the SIL program has shown 
that most of the relevant parameters Influencing Induction times are predicted 
with a reasonable degree of accuracy. The assumption of starting conditions 
where both streams are In chemical equilibrium Is an essential part of this 
effort because (1) the experimental data used for comparison with the program 
results were obtained under conditions where equilibrium species concentrations 
at the combustor inlet Is a reasonable approximation and (2) most of the theo- 
retical approaches for calculating Ignition delays employ the assumption of 
chemical equilibrium as a starting condition. In the next section, the effect 
of non-equilibrium starting conditions Is examined using the SIL program. 




TIME, MICROSECONOS 

FIGURE 7 EFFECT OF WATER VAPOR ON INDUCTION TIMES 


13 









be of significance. The SIL program results of vitiated air and H 2 are shown 
In Figure 10. Non-equilibrium compositions of the vitiated air at the en- 
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entrance Is assumed to have a free-radical composition corresponding to stag- 
nation conditions In the fi manifold. The vitiated air and hydrogen streams 

. . .4 « • tt.A ft X A 3 * 



occur at less than two Inches downstream of the test section entrance. How- 
ever* mixing times are negleeted In this specific analysis* therefore the 
actual delay distance may be somewhat longer than predicted. This figure also 
Includes the delay time calculated assuming chemical equilibrium for both the 

^ S A X » . * . m -A m. «x mm 




spends to a delay distance of approximately 12 inches. 
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FIGURE 10 IGNITION DELAY IOR VITIATED AIR /HYDROGEN 





IV. EXPERIMENTAL PROGRAM 


Sizing of the experimental apparatus was based on the criteria that the 
supersonic Ignition test section should be large, enough to minimize scale 
effects but small enough to reduce run costs. It was also desirable that the 
design conditions simulate a flight regime in which the SCRAM JET vehicle might 
be expected to operate. This regime was chosen so that ignition delay times 
would be relatively long so that they could be measured. 

Design conditions selected for the gas generator were: chamber pressure * 
600 psia, chamber temperature ° 2200°K and flow rate * 9 Ib/sec. The design 
operating conditions for the gas generator and supersonic ignition test section 
are shown schematically in Figure 12. These conditions simulate approximately 
a flight Mach number of 6.5 at an altitude of 80 ,000 ft. , a regime that should 
be of Interest for SCRAMJET application. 

In these experiments, the gas generator pressure and hydrogen pressure 
were selected so that when they expanded to atmospheric pressure the design 
supersonic test section conditions were realized. The combined streams were 
permitted to mix and burn in an unconfined manner as a free jet. The absence 
of a confining wall made it possible to visually observe the flame front and 
determine the ignition delay distance. A solid wall near the reaction zone 
often catalyzes the combustion process by either surface reactions or through 
a temperature Increase due to viscous dissipation. In addition, the flow pat- 
tern (recirculation) becomes more complex. By examining Ignition in a free 
jet, the ignition delay effects should be theoretically predictable providing 
the reaction mechanism is selected properly. 

1. Gas Generato r The components of the gas generator Include an injector, two 
chambers, a turbu 1 ence ri ng and a converging-diverging contoured (parallel flew) 
nozzle. All components of the gas generator except for the Injector were fabri- 
cated out of oxygen free, electrolytic tough pitch copper. A preliminary heat 
transfer analysis indicated than an uncooled copper gas generator system sub- 
jected to the experimental chamber design conditions mentioned previously could 
theoretically operate for run times up to 14 sec. before failure. A maximum 
of 3.5 sec. steady state operation, which entails approximately 6.0 sec. of 
total gas generator run time was necessary for any of the prescribed experi- 
mental runs. 

The two combustion chamber sections are separated by a turbulence ring. 

The purpose of the turbulence ring Is to promote mixing in the hot gases In 
order to insure a relatively uniform exhaust stream. 

The gas generator Injector has three separate elements : The oxidizer in- 
let housing, fuel Injection pintle and a porous natal stainless steel surface 
in the area exposed to the hot combustion gases. The diluent nitrogen is in- 
troduced into the combustion chambers through this porous metal surface and 
thus serves to cool the entire Injector face. The fuel Injection pintle 
serves to introduce the hydrazine Into the combustion chamber and together with 
the nitrogen inlet housing forms an annulus for Injecting the nitrogen tetroxlde 
oxidizer. The pintle contains 10 hydrazine injection ports of 0.0156 in. 
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FIGURE jtt. SCHEMATIC OF EXPERIMENTAL SET-UP 






0 

diameter for the desired flow rate (w f = 1.5 lbs. /sec). 

The gas generator nozzle was designed for uniform, parallel flow across 
the exit plane. Design of the exit contour was accomplished by utilizing a 
computer program which solves the two-dimensional, ax 1 symmetric flow relations 
by the method of characteristics. The program assumes constant values for the 
ratio of specific heats and for the molecular weight of the mixture. This 
assumption yields a good approximation to the contour coordinates calculated 
for real gas conditions. The parallel flow nozzle has a throat diameter of 
1.596 in. and an overall area ratio of 5.38. A photograph of the assembled 
gas generator apparatus is shown in figure 13 with the hydrogen manifold re- 
moved. 

2. Hydrogen Heater The hydrogen heater consists of a large stainless steel 
coll which is heated by convection and radiation from the combustion products 
of three propane burners. It has the capability of supplying conti nuously, 

0.5 Ib/sec of 900° K gaseous hydrogen at an operating pressure of 325 psig. 

3. Supersonic Test Section The exit plane of the gas generator parallel flow 
nozzle together with the hydrogen manifold make up the supersonic test section. 

The hydrogen manifold and nozzle exit section are mated so that an annulus 
is formed between them. The hydrogen gas is expanded to sonic conditions as 
flow through the annulus is choked. The lip separating H 2 and vitiated air 
flow is 0.050 in. thick and was purposely made as small as possible to minimize 
flow recirculation problems in this area, but still maintain some structural 
strength. 

A. EXPERIMENTAL RESULTS 

The experimental program was conducted essentially in two phases. In the 
first phase, the gas generator was developed to a level where it could be oper- 
ated with a high degree of reliability. It was also necessary to demonstrate 
that the oas generator would operate repeatedly at high combustion efficiencies 
(n c > 93%; In order that the gas composition and state properties could be 
reasonably estimated by theoretical techniques. To attempt to achieve high 
combustion efficiency, the gas generator had a large characteristic length 
(L* ^ 120 in.) and a turbulence ring to promote mixing and corrplete combustion 
of the exhaust products. Subsequent to the achievement of satisfactory opera- 
tion of the gas generator, the supersonic combustion ignition delay experiments 
were undertaken. This constituted the second phase of the experimental program. 
Results from each phase of this program are discussed separately in the fol- 
lowing sections. A suiranary of representative performance results from the 
experimental effort, is presented in Table 1. The tabulation presents flow 
rates, chamber pressure and temperature when measured, hydrogen conditions, 
ignition delay data, a characteristic velocity (C*) determined from measured 
values, and a C* efficiency. The theoretical characteristic velocity employed 
to determine the C* efficiency was calculated from an equilibrium thermochemis- 
try computer program. 
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FIGURE 13 GAS GENERATOR ASSEMBLY 



TABLE I 


IGNITION DELAY EXPERIMENTS 


Propellants: N 2 H 4 /N 2 O 4 /N 2 + H 2 


Run No : 

39 

40 

41 

42 

43 

44 

45 

Oxidizer Flow Rate (lb/sec) 

S.62 

5.54 

5.7 

5.7 

5.5 

5.7 

5.5 

Fuel Flow Rate (lb/sec) 

1.58 

1.55 

1.54 

1.54 

1.53 

1.5 

1.5 

N 2 Flow Rate (lb/sec) 

1.98 

2.0 

1.92 

2.04 

1.94 

1*92 

2.0 

Hydrogen Flow Rate (lb/ sec) 

0.28 

0.29 

0.29 

0.21 

0.29 

0.13 

0.08 

f 

G.G. Chamber Conditions 
Chamber Pressure (psia) 

555 

558 

555 

558 

540 

540 

543 

Measured Chamber Temp (°K) 

1680 

1569 

1551 

1534 

1588 

1476 

1476 

Total Flow Rate (lb/sec) 

9.18 

9.09 

9.16 

9/28 

8.97 

9.12 

9.00 

Total Run Time (sec) 

3.3 

3.3 

3.3 

3.3 

3.3 

3.3 

3.3 

C* (Measured) (ft/sec) 

3835 

3893 

3834 

3809 

3835 

3757 

3886 

C* (Theoretical) (ft/sec) 

4105 

4086 

4106 

4067 

4095 

4060 

4052 

C* Efficiency 

0.93 

0.95 

0.93 

0.94 

0.94 

0.93 

0.96 

Measured Ignition Delay (in) 

8 

15 

10 

16 

19 

11 

14 

Theoretical Ignition Delay (in) 

7.6 

8,0 

7.9 

7.97 

8.7 

7.53 

7.19 

H z Measured Manifold Temp (°K) 

750 

645 

710 

691 

590 

681 

705 

Equivalence Ratio 

1.03 

1.09 

1.07 

0.75 

1.12 

0.48 

0.31 
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termlned flame front to be a function of the chemical kinetics of the reacting 
system of gases. There Is also a possibility of Induction of ambient air into 
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psla. Injection conditions for the hydrogen are Mach 1.0 with a static temper- 
ature of 500° K to 800® K (depending on objectives of the run) and a static pres- 
sure of 14,7 psla. The relative velocity ratio for the two streams at Injection 
conditions Is 1 .08 (hydrogen/vltlated-alr) . These experiments, however, had 
relatively large differences In static temperatures for the hydrogen (500°K - 
800®K) and vitiated air system (1000°K) at test section conditions. The rela- 
tive temperature difference for these two streams made It difficult to estab- 
lish an "effective" mean temperature of the gases prior t© Ignition. Therefore, 
an analytical attempt to establish such a temperature was made In order to 
correlate analytical and experimental data. The results of this analytical 
Investigation are discussed In the section entitled, "Discussion of Analytical 








temperature would eliminate this problem. However, If the hydrogen system was 
designed for a static temperature at injection equivalent to the vitiated air 
temperature, the heater would have to be designed to operate at 1200°K. The 
867°K limit on the total temperature for the hydrogen was selected as a com- 
promise to allow the use of relatively inexpensive materials in the construc- 
tion of the hydrogen heater and downstream hardware. 

After system Integrity was confirmed, a subsequent series of experiments 
was initiated to demonstrate the phenomenon of hydrogen ignition in a super- 
sonic (M = 2.9) vitiated air stream. This phenomenon was photographically re- 
corded on black and white 16mm motion picture film taken at 550 frames per 
second. Upon reviewing the pictures, the film velocity (frames per second) was 
sufficient to slow down the ignition process so that an approximate ignition 
delay length could be measured. This measurement is based on the position 
where the flame front (first emitted light) initially appears. A photographic 
history of the Ignition phenomenon is shown in Figure 14. These six consecutive 
frames show the Instant just prior and subsequent to ignition. Frame 1 (upper 
left) depicts the transparent nature of the gas generator products (vitiated 
air). Flow is from left to right. Tn the next frame (lower left) hydrogen 
ignition is assumed to have occurred because of the presence of the visible 
luminous flame front. At this stage in the investigation no deliberate effort 
was made to accurately measure the actual ignition delay length, although it 
was possible to estimate this distance by comparison with known component di- 
mensions in the photograph. A rough comparison of the annular hydrogen injecr 
tor dimensions (10 In. in diameter; with the distance from it to the hydrogen 
flame front, measures an ignition delay length of approximately 10 in. In the 
three runs made during this short exploratory series, the same approximate 
distance was realized. As shown in frames 3, 4, 5, and 6, the hydrogen flame 
front propagates towards the point of Injection and stabilizes approximately 
2 in. downstream. This phenomenon occurred in all of the ignition delay in- 
vestigations* but at slightly longer distances. This stabilized flame front 
can also be considered as a steady state Ignition delay length, but it is felt 
that its position is highly Influenced by the geometry of the hydrogen injector 
lip (recirculation effects). 

In order to completely describe the vitiated air effects on ignition delay, 
as many properties of the gas generator and Its products of combustion (vitiated 
air) as possible must be measured. High temperature 40 percent iridium, 60 per- 
cent rhodium/lrldium thermocouples and water-cooled gas sampling probes were 
fabricated, for use in the measurement of the gas generator chamber temperature 
and the sampling of the vitiated air stream. The purpose of obtaining a vi- 
tiated air gas analysis was to attempt to determine the amount of nitric oxide 
(NO) and nitrogen dioxide (N0 2 ) In the vitiated air, and to compare the actual 
gas composition with that which was theoretically calculated by the Finite Rate 
Reacting Gas (FRRG) program. This specific computer program was discussed in 
the analytical section of this report. Theory and experiments indicate that 
nitric oxide and nitrogen dioxide have a catalyzing effect on the ignition 
process (4). If this Is so, then knowing the mass fraction or existence of 
these two constituents in the vitiated air stream should help justify experi- 
mental results. Also, the presence of nitrogen dioxide in the gas is a measure 
of the Incompleteness of the gas generator combustion process. Nitrogen 
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FIGURE 14 HISTORY OF IGNITION PHENOMENON OF VITIATED AIR AND HYDROGEN 


dioxide is not one of the products of the complete reaction of nitrogen 
tetroxide and hydrazine. 

The analysis of the vitiated air gas sample collected indicated that no 
nitrogen dioxide was present in the air stream. This Implies that the combus- 
tion process in the gas generator must have gone to completion and that the 
inefficiency as manifested in the less than 100 percent C* efficiency was due 
to heat transfer and boundary layer effects. A heat transfer analysis across 
the chamber walls was done and was found to account for at least a 2 percent 
decrease in C* efficiency. 

The heat transfer rate through the gas generator chamber walls was mea- 
sured simultaneously with the chamber temperature and ignition delay lengths. 
This was accomplished by mounting In the chamber wall two chrome! alumel ther- 
mocouples at a depth of 0.813 and 1,563 in. respectively. The distance from 
the deeper thermocouple (1.563 in.) to the hot side of the chamber wall Is 
0.54 in. After each run the slope of the temperature increase for each thermo- 
couple was compared for a quasi-steady temperature difference. The averaged 
difference was 86°K. Knowing the thermal conductivity of oxygen free electro- 
lytic copper and the distance between the thermocouples , it was a relatively 
simple matter of applying Fourier's one-dimensional law of heat conduction for 
a hollow cylinder to get an approximate heat transfer flux rate of 1.54 BTU/ 
in 2 -sec. This efflux of heat was then analytically accounted for in the 
measured C* combustion efficiency of the gas generator. 

To measure the gas generator chamber tenperature two 40 percent iridium, 
60 percent rhodium/iridium thermocouples mounted inside the chamber cavity at 
depths of 0.37 and 1.125 In. respectively from the hot wall, were employed. 

The components of the thermocouple consist of a 0.062 in. inconel shea tit with 
a magnesium oxide layer which Insulates the two 0.010 In. diameter 40 percent 
iridium, 60 percent rhodlum/iridlum wires. A bead joining the two exposed 
thermocouple wires at a distance of 0.25 in. from the sheath material, formed 
the hot junction of the Instrument. The thermocouples held up satisfactorily 
during the course of this investigation. Measurements were taken on all seven 
runs and were relatively consistent as indicated in Table 1. Temperature 
differences between the two thermocouples during the runs when averaged were 
approximately 76*K. The values Of chamber temperature as listed in Table 1 
are the uncorrected readings from the thermocouple positioned closest to the 
hot gas core (1.125 In.). A radiation and conduction thermocouple heat trans- 
fer error analysis was carried out according to the method outlined in 
Reference 12. 

Results of the analysis Indicated an average heat transfer error of 
approximately 5®K and 80®K for radiation and conduction respectively. The 
.primary objectives for measuring the chamber temperature was to verify the 
accuracy of the temperatures computed by the Finite Reaction Rate Gas (FRR6) 
program. Temperatures (^ 1700°K) measured and corrected for thermocouple 
losses did not compare favorably with corrected theoretical temperatures 
(^ 2100°K) calculated by the FRRG program. Reasons for the discrepancy were 
not resolved since the resulting effects of this discrepancy on the investi- 
gation were felt to be minimal. However , due to the many unknowns associated 
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with high temperature thermocouple measurements, the temperature computed by 
the FRRG program were considered to be the more accurate of the two. Therefore, 
the corrected theoretical temperatures were subsequently employed to generate 
an Input temperature which Is used in the Simplified Ignition lag (SIL) program 
to calculate the theoretical Ignition delay distance. The theoretical ignition 
length Is compared with that experimentally measured for correlation. This 
correlation Is discussed In the section "Discussion of Analytical and Experi- 
mental Results." An extensive treatise on this specific portion of the Inves- 
tigation can be found In Reference 13. 

The actual Ignition delay length data were photographically recorded on 
Kodak 4X- Reversal black and white film using a Fastex ISmm camera at a film 
framing rate of 550 frames per second. A photograph of the apparatus for ob- 
serving and measuring the Ignition phenomenon Is presented as Figure 15. The 
graduated measuring scale Is mounted on the side of the test stand for deter- 
mining the Ignition delay length. Since the vitiated air/hydrogen mixture 
burns with a diffuse consistency, the dark background horizontal slot In the 
measuring scale emphasized the Ignition distance and subsequent reaction phe- 
nomena. Each white stripe Is an inch apart. The left side of the horizontal 
slot Is in the exit plane of the hydrogen annular injector and thereby forms 
the datum point from which the exact Ignition delay length Is measured. The 
distance measured from this datum to the point of first light apparition Is 
defined to be the Ignition delay length. Figures 16 and 17 are typical photo- 
graphic data results of the experimental runs. Each figure depicts 14 consec- 
utive frames showing the Instant before and after the onset of Ignition. In 
Figure 16, the fourth frame down shows the point of ignition or a distance of 
10 In. The Inked- In line highlights the line of demarcation between the dark 
background of the horizontal slot In the measuring scale and the vitiated air 
hydrogen combustion flame front. The white elongated marks on the extreme 
outer edge of the film strip are the film velocity timing marks which are gen- 
erated at a frequency of 100 per sec. The film velocity as shown Is approxi- 
mately 550 frames per second. Careful review of the succeeding frames subse- 
quent to the point of ignition, reveals that the flame front propagates towards 
the plane of hydrogen Injection to about 3 in. and then oscillates between 3 
and 7 In. This effect Is believed to be a recirculation phenomenon and a func- 
tion of the nozzle/injector-lip geometry. 

In Run 43 (Figure 17), the temperature of the heated hydrogen (T H * 590K) 

was just slightly above the autoignition temperature. It was demonstrated in 
previous experimental runs that hydrogen Injected at a temperature of 506 °K and 
521 °K would not Ignite In this apparatus . Nevertheless, (Figure 17) Ignition 
was realized but with a relatively long Ignition delay length. Concentrating 
on the third and fourth frames down and In the immediate area to the extreme 
right side of the measuring scale, the onset of Ignition Is Indicated by the 
Illumination of this section. However , the flame front does not sharply appear 
until frame number 6. This Indicates that the point of Ignition Is somewhere 
near the extremity of the measuring scale. By considering the flame front pro- 
pagation speed (approximately 2 In. per frame) and assuming It constant, the 
point of ignition Is calculated to be approximately 19 In. In this particular 
run the flame front stabilized at a distance of 8 in. , which indicates hydro- 
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FIGURE 15 IGNITION DELAY MEASURING APPARATUS 
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FIGURE 16 IGNITION DELAY - RUN 41 
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FIGURE 17 


IGNITION DELAY - RUN 43 









gen tenperature as well as geometry may have a profound Influence on the point 
of stabilization, to attempt was made to resolve this uncertainty. 


Photographic data results of the other four runs are similar In nature to 
these presented. A compilation of the Ignition delay data as a function of hy- 
drogen manifold temperature and equivalence ratio Is presented graphically as 
Figures 18 and 19. A comparison of these results with the analytical findings 
Is discussed In the following section. 
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FIGURE 19 GRAPH IGNITION DELAY VERSUS EQUIVALENCE RATIO 






IV. DISCUSSION OF ANALYTICAL AND EXPERIMENTAL RESULTS 

The ignition delay for a reacting gas mixture is determined primarily by 
the number of free radicals initially present. At no time during the ignition 
and subsequent combustion process does the mass fraction of free radicals be- 
come a large portion of the total gas mixture. Under conditions of relatively 
long ignition delays (low static temperatures), the Ignition process can be 
accelerated significantly by small amounts of free radicals. An extensive 
analysis was undertaken to determine the effect of vitiation on the ignition 
delay characteristics of hydrogen. The analysis showed that the nonequilibrium 
free-radical content of the vitiated air system as It enters the supersonic 
combustion chamber would result in relatively short Ignition delay times. The 
effect is not entirely a function of the vitiated air composition but was shown 
to be related to the nonequilibrium expansion process by which the gas is 
accelerated to test section conditions (recombination freezing). In that sense 
real air exhibits the similar effect on the ignition delay characteristics with 
hydrogen when the air is expanded to the supersonic test section Mach number 
from a high temperature source. 

The nonequilibrium analysis for a test case typical of the experimental 
conditions resulted in a calculated ignition delay length of 1.5 in. for the 
vitiated air. This calculation assumes the hydrogen and vitiated air to be 
instantaneously and completely mixed at a constant mixture temperature of 
1000°K and an equivalence ratio of unity. It also assumes nonequilibrium gas 
compositions for both streams prior to mixing. When equilibrium compositions 
were employed in these calculations, the ignition delay length was computed to 
be 11.5 in. This is approximately a factor of eight Increase In distance over 
the nonequilibrium case. Therefore, should the experimentally observed igni- 
tion delay distance be significantly less than 11.5 In. for mixture tempera- 
tures of 1000°K or less, the effects of the nonequilibrium Inlet conditions 
will have been demonstrated. The fact that ignition occurred in seme cases at 
distances less than 11.5 in., does Indicate the presence of nonequilibrium free 
radicals which have reduced the ignition delay phenomenon. 

An additional theoretical analysis was undertaken to attempt to achieve 
better correlation between experimental and analytical program results. In 
the experimental program, the hydrogen and vitiated air enters the test sec- 
tion area at different static temperatures. This fact complicates the analysis 
of the delay time since a mean mixing temperature for the reacting gases must 
be defined before computing the ignition delay distance. The operational pro- 
cedure for this specific analysis was executed as follows : 

1. In an experimental ignition delay run note the hydrogen manifold tem- 
perature and flow rate, gas generator chamber pressure, propellant flow rates, 
and observed ignition delay distance. 

2. Input the required gas generator properties and parameters into the 
Finite Rate Reacting Gas (FRRG) computer program and calculate velocity , 
specie concentration, static temperature and pressure at the gas generator 
nozzle exit plane, which also serves as the entrance to the supersonic test 
section area. 
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3. Input Into the Simplified Ignition Lag (SIL) computer program the 
specie concentrations computed above, together with the additional hydrogen in- 
jected, the mixture mean temperature, pressure and equivalence ratio. 

4. Determine from the computer- plotted output graph (OH radical) concen- 
tration versus time the point where the slope of constant CM radical orowth 
starts to deviate from straight line behavior (end of Ignition period). 

5. Multiply the time Indicated at this point of deviation by the velocity 
computed by the FRRG to determine the comparable ignition delay length. 

In the initial attempt to compare the experimental and analytical results 
as outlined above, the mean temperature of the mixed stream (hydrogen and viti- 
ated air) was determined by employing the gross over-simplification that the 
mean temperature can be obtained from the weighted heat constant of the com- 
bined streams. Since the flow rate of the vitiated air stream was measured to 
be 9.0 lb/sec in contrast to the flow rate of 0.08 Ib/sec to 0.29 Ib/sec of the 
hydrogen stream, the mean temperature of the mixture, assuming Instantaneous 
and complete mixing, was for all practical purposes equal to that of the 
vitiated air stream (* 980® K) . A typical result of this preliminary effort is 
presented as Figure 20. A comparison of the analytical results with the exper- 
imentally measured ignition lengths were not satisfactory. Experimentally 
there was a factor of approximately 3 difference In the measured ignition delay 
distances for the different hydrogen Inlet conditions. This discrepancy be- 
tween the analytical and experimental results is believed to be entirely an 
effect of the over-simplified assumption of instantaneous and complete mixing 
of the two dissimilar streams i 

A combined mixing and reaction kinetics analysis is beyond the scope of 
the present program, however, a simple rule of thumb used in free jet mixing 
studies states that mixing is usually complete in a distance of approximately 
10 jet diameters from the Injection plane. This distance would be about 45 in. 
In the existing system. The maximum Ignition delay distance observed was 
19 In. *, therefore, it can be safely concluded that the onset of Ignition occurs 
before the hydrogen gas has a chance to diffuse completely into the vitiated 
air core. Hence, the reaction zone of the mixture forms a conical sheath about 
the vitiated air stream. 

In order to improve the correlation of the experimental and analytical 
data, a short analytical investigation of the effects of vi tl ated-ai r/hydrogen 
mixture tenperature on ignition delay as calculated by the Simplified Ignition 
Lag (SIL) computer program was accomplished. Since it was shewn previously 
that the hydrogen does not completely penetrate the vitiated air core before 
ignition, then. It Is only logical to conclude that the mean mixture ignition 
temperature should be a direct function of penetration depth. The Initial 
phase of this specific Investigation employed the assumption that due to the 
limited penetration of the hydrogen into the air core, that ignition occurred 
at a mean temperature equal to the Injection temperature of the hydrogen. The 
static temperature of the hydrogen subsequent to passage through the sonic 
annular injector is the value used In all of the analytical ignition delay 
efforts. Hydrogen temperatures as listed in the tables and graphs are mani- 
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fold or total temperatures. In all cases considered, the OH radical concentra- 
tion never increased. On the basis of this data It can be concluded that 
Ignition never occurred. The experimental observation that Ignition and sus- 
tained combustion did occur In Run 41 and other runs considered In this special 
analysis 9 Invalidates the assumption of mixture temperatures equal to hydrogen 
injection temperatures. In order to approximate the depth of hydrogen penetra- 
tion Into the hot vitiated air core before ignition and. In turn, to be able to 
generate a realistic mean mixture temperature for analytical purposes, the 
GASL's "Finite Rate Supersonic Combustion Mixing Analysis" (14) computer pro- 
gram was employed. This program describes the turbulent mixing of axisymrrstric 
hydrogen-air jets inside a duct. Even though the program was written for 
ducted flow, the Initial portion of the output Is applicable to the mixing and 
Ignition processes. This is only true when the processes are confined to a 
region located a considerable distance upstream from the point where the reac- 
tion zone would theoretically contact the duct wall. 

The approximation of a mean mixture temperature was obtained by using the 
calculated penetration depth at which the equivalence ratio output of unity 
occured. This value was located at a grid point positioned 0.414 in. deep 
into the periphery of the vitiated air core at an axial distance of 1.2 In. 
from the plane of hydrogen Injection. It was therefore concluded that in the 
distance from the hydrogen Injection plane out to this station (1.2 In.) the 
hydrogen gas diffused Into the vitiated air core to a depth of 0.414 in. 

The volume percentage of the vitiated air core penetrated by the hydrogen gas 
in this distance amounts to approximately 40 percent. Therefore, the mean 
temperature of the mixture was arrived at by using only 40 percent of the 
total vitiated air flow In a v1tiated-a!r/hydrog@n weighted mean temperature 
calculation. This effort reduced the mean mixture temperature to a more 
realistic value and resulted In a better correlation between analytical and 
experimental data. The static temperature of the vitiated air core at the 
hydrogen Injection plane was calculated to be 850°K. This value was arrived 
at by subtracting the temperature drop (AT *» 133®K) due to heat transfer 
effects In the conbustlon chamber and along the nozzle walls, from the theo- 
retical static temperature (983®K) calculated by the FRRG program. In the 
heat transfer analysis. It was found that the chamber temperature decreasedgby 
67°K for a quasi-static heat transfer rate through the walls of 1.55 BTU/ln - 
sec. It was therefore assumed that the additional temperature decrease through 
the nozzle section should be at least equal to that lost In the combustion 
chamber. With a modified vitiated air static temperature of 850°K versus the 
original theoretical static temperature of 983°K 8 It Is felt that the species 
concentrations as calculated by the FRRG program at the original temperature 
(983°K) would not vary excessively due to this temperature difference. Hence, 
the species concentration as calculated originally at 983°K was input into the 
SIL program along with the mean temperature generated by using the modified 
temperature (880°K) . A typical result representative of this analysis is 
shown in Figure 21 which is manifested In the longer Ignition delay time. It 
was also observed that the analytical and experimental Ignition delay lengths 
agreed much better than In the previous comparison where It was assumed that 
the hydrogen and vitiated air streams mixed Instantaneously and completely 
prior to Ignition. This comparison would seem to Indicate that the effect of 
different static temperatures for the hydrogen and vitiated air streams In 
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addition to mixing rates , have a profound effect on any analytical ignition 
delay model. 

Figures 22 and 23 are graphical plots comparing the experimental and ana- 
lytical data as a function of hydrogen manifold temperature and equivalence 
ratio. As can be observed In Figure 22, the differences In delay length diverges 
as the hydrogen temperatures decrease. This effect may be caused by the effect 
of temperature on the mixing or diffusion rate in the two streams. It is shown 
in Reference 15 that the mass diffuslvity In gaseous mixtures does increase with 
temperature. Another noteworthy observation Is the trend of decreasing Ignition 
length with increasing hydrogen temperature as manifested In both the analytical 
and experimental data results. This fact tends to reconfirm the validity of the 
analytical model in that Ignition delay length decreases with increasing tempera- 
tures. In Figure 23 the experimental data deviate from the analytical theory of 
relative Independence of Ignition delay distances on equivalence ratio. The 
observed fact that the experimental data points do straddle a straight line 
when drawn equidistant from each point, provides some encouragement. An expla- 
nation as to why there Is such scatter cannot be given at this time. Analy- 
tical! y, the data points do Indicate a relatively constant Ignition delay dis- 
tance Irregardless of equivalence ratio variance. 
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FIGURE 23 GRAPHICAL COMPARISON OF ANALYTICAL AND EXPERIMENTAL IGNITION DELAY 
RESULTS AS FUNCTION OF EQUIVALENCE RATIO 




V. CONCLUSIONS 


Detailed analytical and experimental measurements have been made In a 
mixing and reacting vitlated-alr/hydrogen supersonic stream for the purpose of 
determining the effects of vitiation products (free radicals, third bodies, 
etc.) on the Ignition delay phenomenon. On the basis of the experimental data 
and their subsequent correlation with the analytical model, it is concluded 
that: 

1. The presence of water vapor In a combustor entrance flow decreases the 
ignition delay or induction time for the temperature regime Investigated. 

2. The analytical ignition delay model was not independent of the mixing 
process. Results predicted by the analytical model compared favorably with the 
experimental data when effects of mixing were considered. Therefore in practi- 
cal systems, reaction and mixing are coupled and should be considered simul- 
taneously. 

3. Nonequilibrium free-radical content of the Incoming stream accelerates 
the reaction process. This causes Ignition to occur much sooner than would be 
expected if the incoming stream were near equilibrium as far as free-radical 
concentrations are concerned as occurs when atmospheric air is diffused to the 
supersonic combustor Inlet conditions. This Indicates that the use of vitiated 
air to simulate the ignition phenomenon for a SCRAM JET system will lead to data 
which differs from flight conditions. These conclusions apply in the low tem- 
perature supersonic regime where reaction times are rate controlling. 

4. In the temperature regime Investigated (* 850°K) , the experimentally 
measured content (4 percent) of nitric oxide (NO) in the vitiated air had little 
effect on the ignition delay length. The catalytic effect of NO on ignition 
was not observed. 
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